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Abstract 

The properties of the high energy pions observed in heavy ion collisions, in particular 
in the system Au on An at 1 GeV/nucleon are investigated. The reaction dynamics is 
described within the Quantum Molecular Dynamics (QMD) approach. It is shown 
that high energy pions freeze out early and originate from the hot, compressed 
matter. A*-resonances are found to give an importnat contribution toward the high 
energy tail of the pion. Further the role of in-medium effects in the description of 
charged pion yield and spectra is investigated using a microscopic potential derived 
from the Brueckner G-matrix which is obtained with the Reid soft-core potential. 
It is seen that the high energy part of the spectra is relatively more suppressed 
due to in-medium effects as compared to the low energy part. A comparision to 
experiments further demonstrates that the present calculations describe reasonably 
well the neutral (TAPS) and charged (FOPI) pion spectra. The observed energy 
dependence of the ratio, i.e. deviations from the isobar model prediction, 

is due to Coulomb effects and again indicate that high energy pions probe the hot 
and dense phase of the reaction. These findings are confirmed independently by a 
simple phase space analysis. 
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1 Introduction 


Heavy ion collisions play a vital role in the quest for the nuclear equation 
of state, since a hot and dense nuclear matter phase is created, though tran¬ 
siently, in heavy ion collisions. In a typical intermediate energy reaction, for 
e.g. Au-|-Au at 1 GeV/A, matter gets compressed to densities of about two to 
three times the normal nuclear matter density po and the maximum temper¬ 
ature attained in this central zone is supposed to reach values of few tens to 
about 100 MeV [1-3]. 
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The investigation of particle production is a well established approach to probe 
the nuclear matter under these extreme conditions [4-6]. At incident energies 
around 1 GeV/A the pion production is the dominant channel by which the 
hot, compressed matter de-excites. Pion spectra at these relativistic energies 
have been measured at GSI (FOPI [7], KaoS [8,9], TAPS [10,11]) and LBL 
[4,5]. Unlike at lower energies the spectra show non-thermal slopes, i.e. one 
needs a superposition of two Maxwell-Boltzmann distributions to describe the 
data [5,7,9,11-13]. The onset of the second slope roughly coincides with the 
threshold for the A—resonance. Therefore, the resonance production and its 
subsequent decay into pions is considered to be important for an understanding 
of the excited nuclear matter. For this reason, in addition to the A—channel we 
include one-pion and two-pion decays of A*-resonances and the respective one- 
pion reabsorption process. The role of these resonances in the production of the 
high energy pions is analysed in this work. Furthermore, recent measurements 
indicate that the high energy pions freeze out early, particularly when emitted 
perpendicular to the reaction plane [8,14]. Therefore, these pions most likely 
probe the early phase of the collision. Hence, a study of subthreshold pions is 
of current interest and is being actively pursued by both experimentalists and 
theoreticians. 

As compared to the nature of subthreshold pions, the properties of the low 
energy pions, which are interpreted to be originating from A—decay, are rela¬ 
tively well understood within transport model approaches. However, one hnds 
that the transport models such as QMD and BUU, in general, overestimate 
the charged pion yields as well as the low energy part of the spectra [7,15]. In 
this context, one has to keep in mind that the nuclear mean held, in particular, 
the respective momentum dependence of the interaction has large effects on 
various dynamical observables [16-18]. In most of the previous studies on pion 
production a nuclear mean held derived from ehective phenomenological inter¬ 
actions, such as the Skyrme forces, have been applied [3,7,15,19,20]. However, 
a more realistic description of baryon dynamics in the hot and compressed 
phase of nuclear matter is possible within the framework of Brueckner theory 
where the nuclear mean held is derived microscopically from the G-matrix ob¬ 
tained by solving either the Bethe-Goldstone or the Bethe-Salpeter equation. 
With these realistic forces one hnds in general an enhancement of the repulsive 
part of the nuclear interaction which originates from the strong momentum 
dependence of the G-matrix [16,17,21,22]. Since observables such as nuclear 
how are found to be sensitive to the use of realistic forces [17,23,24], it is a 
worthwhile task to analyse the role of in-medium ehects in the description 
of pionic observables at SIS energies. Therefore, in the present work, we also 
investigate the ehect of the Brueckner G-matrix potential on pion multiplici¬ 
ties and spectra and then compare to the results obtained with the standard 
Skyrme interaction. 

Another interesting aspect of the pion production observed in heavy ion col- 
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lisions is its charge dependence [25-29]. In the case of heavy systems such 
as An one finds [7,9] for the total abundancies as a general law > 

The two principal effects contributing to the observed charge 
dependence are the isopin dependence of the pion production cross sections 
and the Coulomb force. Since most of the pions observed at bombarding en¬ 
ergies around 1 GeV/A originate from the decay of A—resonances, one can 
determine the isospin effect in a simple isobar model [6,30,31], which depends 
only on the Z/A ratio of the colliding systems and is independent of the pion 
energy. However, in Refs. [7,9] deviation from the behaviour predicted by the 
isobar model have been observed. The observed energy dependence of the 
7r“/7r+ ratio is then essentially due to the other factor, i.e. the Coulomb effect 
which modifies the pion momentum and the available phase space [28]. As 
mentioned in Ref. [32] one can moreover use the energy dependence of this 
ratio to estimate the size of the pion emitting source and thereby the density 
of the participant zone probed by the observed pions. This provides us with an 
alternative approach to investigate the nature of the subthreshold pions. To 
study this aspect, we take into account the Coulomb interaction between the 
baryons and the charged pions. This contributes toward the real part of the 
pion optical potential. Another contribution arises from the effect of nuclear 
medium on the pion dispersion relation, i.e. from collective AN~^ and NN~^ 
excitations. However, it is shown in a recent study [33] that only weak correc¬ 
tions, i.e. a small in-medium pion-nucleus potential, are required in order to 
obtain a realistic description of pionic observables. In addition, also in other 
works [34,35] it was found that the high energy part of the pion spectrum is not 
much affected by these medium corrections. A more sophisticated treatment 
of collective pionic excitations within the transport approach was proposed 
in Ref. [36] which was subsequently applied to heavy ion collisions in Ref. [37]. 
However, as discussed in Ref. [38], it is rather questionable if the simple AN~^ 
( and NN~^) model yields a readable pion dispersion relation. Further, the 
results of Ref. [37] indicate that the influence of collective excitations on pio¬ 
nic observables is small and tends to enhance the pion yields with respect to 
the standard approaches and the experimental findings. Hence, we omit such 
colective medium effects on the charge dependence of pion production in the 
present work. 


Here we mainly focus upon the aforesaid aspects and investigate them using 
the Quantum Molecular Dynamics (QMD) model [39]. The essentials of the 
QMD model is briefly introduced in section 2. Results obtained are discussed 
in detail in section 3. Finally, we summarize our hndings in section 4. 
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2 QMD model 


Quantum molecular dynamics is a semiclassical transport model which ac¬ 
counts for relevant quantum aspects like the Fermi motion of the nucleons, 
stochastic scattering processes including Pauli blocking in the hnal states, 
the creation and reabsorption of resonances and the particle production. A 
detailed description of the QMD approach can be found in Refs. [17,39]. 

The time evolution of the individual nucleons is thereby governed by the clas¬ 
sical equations of motion 


dpi _ _d^ ^ 

dt dqi ’ dt dpi 

with the classical A^—body Hamiltonian H 

H = T. + 5 E {u.s + + C'F"') ■ (2) 

i DJ 


The Hamiltonian, Eq. (2), contains mutual two- (and three-) body potential 
interactions which are hnally determined as classical expectation values from 
local Skyrme forces Uij supplemented by a phenomenological momentum de¬ 
pendence and an effective Coulomb interaction 


Uij = a + f3 + 3ln^ (e\pi - -h l) —, 

\Po J \Po J ^ ^ Po 

Wi - 


TJCOUl _ / _ \ ^ 

~\a) %-q, 


-erf 


where pij is a two-body interaction density dehned as 

1 


Pij 


(47rL) 




( 3 ) 

( 4 ) 

( 5 ) 


and erf is the error function. The parameters a,/?,7,5, e of the Skyrme in¬ 
teraction, Eq. (3), are determined in order to reproduce simultanously the 
saturation density (po = 0.16/m“^) and the binding energy {Eb = —16 MeV) 
for normal nuclear matter for a given incompressibility as well as the 
correct momentum dependence of the real part of the nucleon-nucleus optical 
potential [17,40]. In the present calculations we use a soft equation of state 
with Koo = 200 MeV including the momentum dependent interaction (SMD), 
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Eq. (3). The Yukawa-type potential in Eq. (2) mainly serves to improve 
the surface properties and the stability of the initialized nuclei. 

Nucleons while propagating shall collide stochastically, if the distance be¬ 
tween the centroids of the two Gaussian wavepackets is less than dmin = 
\/o'tot(v^)/7r. This collision process is implemented using Monte Carlo method 
and the effect of Pauli exclusion principle is duly taken into account. For 
the inelastic nucleon-nucleon channels we include the A(1232) as well as the 
TV* (1440) resonance. In the intermediate energy range the resonance produc¬ 
tion is dominated by the A, however, the N* yet gives non-negligible contribu¬ 
tions to the high energetic pion yield [41]. The resonances as well as the pions 
originating from their decay are explicitly treated, i.e. in a non-perturbative 
way and all relevant channels are taken into account. In particular we include 
the resonance production and rescattering by inelastic NN collisions, the one- 
pion decay of A and N* and the two-pion decay of the N* and one-pion 
reabsorption processes. ( For details see Ref. [33].) 

Pions thus produced are propagated according to the same equations of motion 
given in Eq.( 1). The Hamiltonian corresponding to pions is given as 

/^. = E(Al + ™5 + r.°°"')- (6) 

i=l ^ ^ 


The Coulomb interaction between the baryons the charged pions is given 
as 


Nh 




Coul 


= E 


€i€j 


=i \ Qj - Qi 


(7) 


where Ni, and are respectively the total number of baryons including 
charged resonances and pions at a given time t. Thus, pion propagation at 
the intermediate and hnal stages is guided essentially by the Coulomb effect. 


3 Results and Discussions 


3.1 High energy pions and role of N* resonances 


Since the observation that high energy pions freeze out early when emitted 
perpendicular to the reaction plane, the focus has been to establish, both 
experimentally and theoretically that these high energy pions probe the hot 
and dense phase of the collision. Hence, we analyse the properties of high 


5 



energy pions and the effect of heavier baryon resonances, such as the iV*(1440), 
in this regard. Here we adopt the phenomenological Skyrme interaction SMD. 

For this purpose we calculated the total pion multiplicity in Au+Au collision 
at 1 GeV/A corresponding to the minimum bias condition. According to Ref. 
[7], we choose 6niax=ll fm. Results obtained are shown in Fig.l as a function 
of pion transverse momentum pt and pion freeze out time tfreeze, he. the time 
after which pions do no more interact via inelastic scattering processes. It 
can be seen that most of the pions are in the region pt < 200 MeV/c and 
tfreeze ~ 30 fm/c. lu Order to project out the high energy pions, we normalised 
the pion multiplicity M{pt) obtained for each value of pt by its maximum value 
Afmax(pt)- The normalized multiplicity is shown in Fig.2, where it is seen that 
the peak value shifts toward early times as the value of pt is increased. For 
example, pions with pt ~ 200 MeV/c freeze out at t ~ 30 fm/c; whereas 
for Pt ~ 500 MeV/c one finds tfreeze ~ 20 fm/c. This directly shows that 
high—Pi pions originate from early times. In order to relate these pions with the 
compressed early phase, we calculated the total baryon density ps in a sphere 
of 2 fm radius around the collision centre with impact parameter b=l fm. The 
time evolution of the ps/Po is shown in Fig.3, where the po = 0T6 is the 

normal nuclear matter density. The maximum density is achieved at about 25 
fm/c and is about 2po. Over the time span 10 < t < 30 fm/c the density in the 
collision center is at least po, and this time span may then be identified with 
the hot, compressed participant phase. Thus, from this analysis we find similar 
to a previous study using hard equation of state and only A—resonances [20] 
that high energy pions with ffreeze ~ 20 fm/c pertain to the early phase of the 
collision, when the matter is compressed to about 2po- 

For the sake of completeness we show in Fig.4 for two impact parameters, 
b = 0 fm and 6 = 6 fm the collision configuration, i.e. the evolution of 
the average distance AR = Rp — Rp between the centres of the colliding 
nuclei where Rp and Rp represent the average position of the centres of the 
projectile and target nuclei, respectively. The collision centre is at the zero of 
the Y-axis. The curves representing Rp are in the positive Y-axis plane and 
those representing Rp are in the negative Y-axis plane. The two dotted lines 
correspond to the sharp density radii of the two colliding nuclei. Thereby, it 
can be seen that at about 4 fm/c when the curves corresponding to Rp and 
Rp intersect the dotted lines the nuclei come in contact, i.e. this represents the 
touching configuration. Further, for very central collisions the time span upto 
~ 25 fm/c pertains to the compression stage, and beyond which the expansion 
process begins. Hence, the maximum density or compression is attained at 
about 25 fm/c, which is consistent with the results shown in Fig.3. It may 
also be noted that for finite values of impact parameter b, one has at the time 
of maximum compression AR = b. 

To illustrate the contribution of high energy pions from N* resonances, we have 
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compared in Fig.5 the total pions originating from A— resonance to those from 
N*. The pion multiplicity linearly weighted by impact parameter b is shown as 
a function of b. It is seen that the contribution from iV*-resonances irrespective 
of the nature of the collision, i.e. central, semi-central or peripheral, is about 
5%. Though their multiplicity is small, they play an important role as most 
of them originate from early times. This aspect is illustrated in Fig.6 where 
we have plotted as a function of pion freeze out time tfreeze- The total 

number of pions from A and N* are labelled as and M^* respectively, 
with ■ Interestingly, the maximum contribution from N* 

(about 8%) occurs at tfreeze ~ 20 fm/c. It maybe recalled here that pions with 
this value of tfreeze pertaiu to the dense phase with pB > 2po- Therefore, the 
majority of the N* pions originate from early times and hence are most likely 
high energetic ones. 

To demonstrate this aspect we calculated the neutral pion spectrum at midra¬ 
pidity under the minimum bias condition. ( We have chosen tt^ so that the 
high energy tail of the spectrum is not distorted by Coulomb effects.) Results 
obtained are shown in Fig. 7 alongwith the 7r° spectrum calculated only with 
N* pions. In the same hgure, we have also shown the TAPS data for compar- 
ision. Due to a recent experimental analysis the data given in [11] have to be 
renormalized by factor 0.6 [15]. It is clearly seen that the N* contribution to 
the high energy tail is by a factor of about two more than its contribution to¬ 
ward low energy pions. This is in agreement with a recent study [15] based on 
the BUU model where it was shown, though for a different system at a higher 
energy, that contributions from baryonic resonances like A* (1440,1520,1535) 
and A(1600) are significant over the high energy tail of the pion spectra. More¬ 
over, the 7r°(A*) spectra in Fig. 7 portrays two distinct slopes and the slope 
parameter T extracted from the high energy part is relatively higher than the 
one obtained from the high energy part of the total tt® spectrum. This again 
establishes that N* pions pertain to earlier times as compared to those from 
A. Therefore, it may be said that pions from N* need to be included for a 
consistent description of the high energy tail. 


3.2 In-medium effects on charged pions 


As discussed in Introduction, we now investigate the effect of microscopically 
determined nuclear mean field on obsevables such as pion yield and spectra. 
(All data shown in this section are taken from Ref. [7]. ) 

For this purpose, the G-matrix is obtained by solving the Bethe-Goldstone 
equation 


Q 

G = !/ + !/-^- G 

u - Ho-\-ir] 


( 8 ) 


7 



where Q is the Pauli operator which forbids nucleons from scattering into 
occupied states. The bare NN interaction V is taken to be the Reid soft-core 
interaction and Hq is the single particle energy [42,43]. Thus, the G-matrix 
accounts for the most relevant medium effects in nuclear matter. The mean 
held potential Ggmat is derived from the G-matrix as [16,17] 

Ggmat(p,:^ = 4 / I G{p) I p,^) . (9) 


For the purpose of a simple application within the QMD model the potential 
is parametrized in a Skyrme-like form 

p /9 ^3 o o /9 ^6 

UcMATip^p) = ai{—) + a2{—) + ajn^la^p^ + !]{—) . ( 10 ) 

Po Po Po 


The coefficients a's were determined from a ht of Eq. (10) to the micro¬ 
scopically calculated values of Ggmat, Eq-(9), within the density range of 
0 < p < 2po and for momenta p < 500 MeV/c. The momentum dependence 
of so-determined potential at p = po is shown in Fig.8 alongwith the one 
obtained from the phenomenological Skyrme force using Eq. (3). (Galcula- 
tions done with the Skyrme force and the G-matrix potential are hereafter 
referred to as SMD and GMAT, respectively.) For comparison, we have also 
shown the potential as obtained by Li and Machleidt [21] using the Dirac- 
Brueckner approach. It can be seen that the nuclear mean held derived from 
the non-relativistic G-matrix is slightly more attractive at low momenta and 
less repulsive at high momenta compared to the one obtained from the rela¬ 
tivistic G-matrix. However, the general behaviour is quite similar. The promi¬ 
nent feature of the realistic G-matrix interaction is the signihcantly stronger 
momentum dependence as compared to the Skyrme force. In a recent study 
[44] on the ehect of realistic forces on collective nucleon how it was found that 
the relativistic (Dirac-Brueckner) and non-relativistic (Brueckner) approaches 
yield quite similar results when applied to heavy ion collisions. However, the 
diherences to phenomenological interactions (Skyrme or Walecka model) are 
signihcant. Thus, it is of interest to investigate the sensitivity of pionic ob¬ 
servables to realistic forces. In order to be as consistent as possible, we have 
also considered the ehect of in-medium cross sections [16,17] in the elastic NN 
channel for center-of-mass momentum p < 500 MeV/c, which are derived from 
the same microscopic G-matrix potential using the two-Fermi-sphere aprroach. 
However, the ehect of the in-medium cross section on the pion yield has been 
found to ba small compared to that of the mean held and is therefore not 
discussed in detail in the following. 

Results obtained for the 7r~^ multiplicity is plotted in Fig.9 as a function of 
the participating nucleons Apart- R is seen that the theoretical calculation 



overestimate the data of the Au system by about 50% at small impact param¬ 
eters. In this context we want to mention that the pion multiplicities obtained 
by other groups for central Au + Au reaction at 1 GeV/A, using QMD [3,7] 
or BUU [19,37] models, differ only within about 10% with our results. Thus, 
this overestimation at small impact parameters seems to be a general feature 
of present transport calculations. However, the theoretical prediction is more 
in agreement with the experimentally derived Harris systematics [45]. It may 
also be noted that the Harris systematics while it overpredicts as compared 
to the FOPI data in the case of heavy system such as Au is in good agree¬ 
ment with the data for smaller system like Ni [46]. We have also shown in 
Fig.10, the rapidity distribntion obtained nnder minimnm bias condition. In 
both SMD and GMAT calculations one hnds at midrapidity a overestimation 
of abont 30%, as compared to the data. Similar featnres are also observed in 
the case of n~. 

To see the effect npon momentnm distribntion, we calculated the charged pion 
spectra. Results obtained with both SMD and GMAT potentials are then com¬ 
pared in Figs.11 and 12 to the corresponding FOPI data. The overall agree¬ 
ment with the observed spectra is reasonably good. In particular, one hnds 
that with the inclusion of in-medium effects, there is a slight reduction over the 
low energy part thus coming closer to the data. Bnt, a larger effect is observed 
over the high energy part where the snbthreshold pions get qnite snppressed 
indicating the importance of resonances heavier than A and N*. Thus, one 
hnds that, as compared to the high energy part, the low energy part is not too 
mnch ahected by the strong momentnm dependence of GMAT potential which 
may be dne to the fact that low energy pions are freqnently reabsorbed. For a 
more complete description of pion yields and spectra, one probably needs to 
inclnde the in-medinm ehects in the inelastic NN channel as well. In addition, 
a possible reason for the overestimation over the low energy part can be the 
life-time of the resonances involved. ( Some indication towards this can be 
inferred from Fig.l, where one hnds that even at 60 fm/c, resonances are still 
decaying into pions thongh their mnltiplicity is qnite small. ) E.g. in Ref. [47] 
an alternative treatment of the resonance lifetime has been proposed, which is 
based on more qnantnm mechanical aspects. However, it is still nnclear at the 
moment how snch a description should be applied consistently in transport 
simnlations. 


3.3 TT /7r+ ratio and Coulomb effect 


We had demonstrated theoretically that high energy pions originate from the 
early phase of the collision. Here, using the observed energy dependence of 
7r“/7r+ ratio, we re-iterate that high energy pions do probe the hot and dense 
matter. 
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It is known that the two principal effects contributing toward the observed 
charged dependence of pion production, i.e. are 

the isospin dependence of the pion production cross sections and the Coulomb 
force. As most of the pions observed are from A—resonance, one can determine 
the effect of isospin. Using the isobar model [6] one hnds for pion originating 
from A-decay the general relation 


a{7r )/(T(7r°) = 


5A2 + NZ 

AT2 I AAT'y 




5^2 + NZ 


/ -I -I \ 


For the Au+Au system we obtain the ratios tt” : 7r° : tt’*' = 1.37 : 1 : 0.71 for 
the repective multiplicities. Therefore, one has = 1.95. For pions from 

N* the isobar model yields = 1.70. Since these relations are globally 

valid, i.e. independent of the pion energy, the observed energy dependence 
of the 7r“/7r+ ratio must be due to the Coulomb potential energy Vq [20,19]. 
This is illustrated in Fig. 13, where the ratio of the calculated differential cross 
sections with ge = d2(T/(pEdp dhl) are compared with the 

observed one [7]. One hnds that the overall behaviour of the measured ratios 
is well reproduced. In particular, over the low energy part both SMD and 
GMAT curves are quite similar. However, at higher energies T^.m. > 300 MeV, 
he. Pc.m. > 400 MeV, the GMAT results show stronger huctuation as compared 
to the SMD case, which may be due to the suppressed cross sections. Such 
strong huctuations over high energy part are also apparent in the experimental 
data. Nevertheless, the trend in Fig. 13 indicates that the results approach a 
constant value as the pion momentum Pc.m. increases which is in agreement 
with the experimental Endings. 


The general features of the efiect of the Coulomb force on the can be 

understood using a simple analytical model as described in Ref. [28]. Let po 
and p be the initial and final momenta of a charged particle. Then the charged 
particle cross section a{p) = d^a/d^p can be related to the neutral particle 
cross section ao(p) = d^Uo/ d^p through the change of 

(^{p) = Mpo{p)) I ^ I • (12) 


The two effects due to Coulomb force are the momentum shift 6p = p — po 
and the Coulomb phase space factor | d^po/d^p |. Taking relativistic quantum 
corrections up to first order in Za and in the limit p —> oo, one has ±(5p ~ 
Vc/Pp and | d^po/d^p |~ I ±Vc/p, where Pp = p/ojp and ujp = ^p^ + m\. 
Omitting the phase space factor in Eq. (12) in the limit p —> oo the 7r~/7r~^ 
ratio is then given as 

( 13 ) 
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where uj‘^= p\ + ml and p± = pT ^P- 


It is found experimentally that the high energy slopes of 7r“ and tt’*' are very 
similar and hence the slope parameter T is the same for both positive and 
negative charged pions. Moreover, one hnds Eq. (13) is independent of the mo¬ 
mentum p and the ratio 'n~ is in good approximation given by a constant 
determined by the Coulomb energy shift, i.e. by Vq- This is consistent with 
the experimental observation where one hnds that the ratio approaches 

a constant value of about 0.93 in the high momentum limit. Therefore, the 
above approximate equation can be used in the high energy limit to estimate 
the Coulomb energy Vc- In addition we take the value of C_/C+ according 
to Eq. 11 to be 1.95 thus accounting for the isospin dependence of the cross 
sections. It may be noted that a very similar value (C'_/C'+ = 1.94 ± 0.05) is 
obtained experimentally in the zero momentum limit [32]. Taking the slope 
parameter T to be 94.6 MeV [7] and the constant value of ti ~ratio to be 
0.93 we obtain the Coulomb energy Vq to be about 35 MeV. The radius R of 
the pion emitting source is then estimated to be 4.5 fm with the assumption 
of a static uniform Coulomb held Vc ~ Ze^/R where the total number of 
participating charges Z is taken to be 110 as estimated in Ref. [32]. We also 
estimated the density p ~ 3A/(47ri?^) of the participant zone to be about 4po 
where the total number of participants A ~ 197 {ZjlA). It may be noted that 
the estimation of p is about two times the value obtained from the actual 
simulation, see also Fig. 3. If we consider a more realistic value p = 2po and 
using the fact that Vc — 35 MeV we obtain R = 6.7. This then gives an idea 
about the order of magnitude of uncertainity present in the estimation of the 
pion source size R and the effective number of participating charges Z. From 
this simple analysis based on the observed energy dependence of the n~/n~^ 
ratio, we hnd in support of our theoretical investigation that the high energy 
pions probe the hot, compressed phase of the collision. However, the empirical 
determination of the freeze-out radius from measured ratios [32] should 

be performed carefully. In this context, it may be noted that a recent analysis 
[29] of Ti ~ratios measured in Au + Au collision at 11 GeV/A and in 
Pb + Pb at 158 GeV/A shows valley like structures in the minimum plane 
rendering the determination of R and Z from the data difficult. Moreover, 
the scenario for the origin of high energy pions is probably not that of an 
equilibrated and thermalized source but rather that of two counterstreaming 
nuclear matter currents with high relative velocities and thus, e.g., Lorentz 
forces should be taken into account in such an analysis. 


4 Summary 


In this work we investigated the production of pions in Au-|-Au collisions at 
1 GeV/A using the QMD model. Focussing upon the high energy pions, we 
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showed that these subthreshold pious freeze out at early times when the matter 
is compressed to more than two times the saturation density. Hence, they probe 
the hot and dense phase of the nuclear reaction. One of the processes which 
can produce such energetic pions is the decay of baryonic resonances like the 
N*. A detailed analysis shows that eventhough the multiplicity of pions from 
the decay of N* is small compared to those from A—decay, i.e. about 5% of the 
total yield, they play an important role as most of them pertain to the early 
phase of the reaction. Further, it is found that N* contribution toward high 
energy pions is by a factor of about two more than its contribution toward 
low energy pions. Therefore, pions from higher resonances like N* need to be 
included for a consistent description of the high energy tail of the pion spectra. 

Moreover, due to the importance of momentum dependence in the study of 
heavy ion collisions, we have also studied the role of in-medium effects by per¬ 
forming calculations with a microscopic nuclear mean held derived from the 
Brueckner G-matrix using the Reid soft-core interaction. It is found that cal¬ 
culated pion multiplicities, though they agree reasonably well with the Harris 
systematics overestimate the FOPI data at small impact parameters by about 
50%. The use of realistic G-matrix interaction improves this situation, how¬ 
ever, not to a sufficiently high extent. We then compared the charged pion 
spectra calculated with G-matrix and Skyrme potentials to the respective 
FOPI data. It is found that the high energy part gets relatively more sup¬ 
pressed by in-medium effects, as compared to the low energy part. Therefore, 
for a complete description of pion yield and spectra, one probably needs to 
consistently include in-medium effects both in the elastic and inelastic NN 
channels. In addition, an improved treatment of lifetime of resonances may 
also be warranted. 

Finally, we investigated the isospin dependence of the pion production. The 
present calculations are able to reproduce the overall behavior of the mea¬ 
sured ratios. Deviations from the simple isobar model originate from 

the Goulomb force. These deviations found for high energy pions indicate the 
influence of the Goulomb held on pions stemming from the high density phase. 
These hndings are in qualitative agreement with a simple analysis based only 
on the observed 7i~/7i~^ ratios and the slopes of the respective spectra. 
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Fig. 2. Normalized pion multiplicity, i.e. the multiplicity calculated for each value 
of pt is normalized by its maximum value, obtained in as the same reaction in Fig. 
1 is shown as a function of pion transverse momentum pt and pion freeze out time 

^freeze • 
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Fig. 3. Time evolution of the baryon density /jb in a central ( b=l fm) An on Au 
reaction at 1 GeV/A. pB is calculated with a Skyrme force (SMD) in a sphere of 2 
fm radius around the collision center. 
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Fig. 4. Time evolution of the average distance AR = Rp—Rx between the centres of 
the two colliding nuclei {Au + Au) is shown for two impact parameters. Rp and Rp 
represent the average position of the centre of the projectile and target, respectively. 
Curves corresponding to Rp are in the positive Y-axis plane, and those for Rp are 
in the negative Y-axis plane. Zero of the ordinate corresponds to the collision centre. 
Dotted lines represent the sharp density radii of the two nuclei. 
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Fig. 5. The total multiplicity of pious calculated with a Skyrme force (SMD) orig¬ 
inating from A^*-resonances is compared to those from A-resonances. The pion 
multiplicity linearly weighted by the impact parameter b is shown as a function 
of b. 
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Fig. 7. Neutral pion pt spectrum calculated with a Skyrme force (SMD) at midra¬ 
pidity for a Au on Au reaction at 1 GeV/A under minimum bias condition is com¬ 
pared with the corresponding tt® spectrum of pions originating from the decay of 
N*-resonances. For comparison, TAPS data given in Ref. [11] is also shown. Due 
to a recent experimental analysis the data are renormalized by factor 0.6. 
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Fig. 8. The momentum dependence of nuclear single particle potential derived from 
the Brueckner G-matrix (GMAT) using the realistic Reid soft-core NN interaction is 
compared to the Skyrme interaction (SMD) at a density p = po- For comparison, also 
the potential obtained in Ref. [21] using the relativistic Dirac-Brueckner approach 
(DBHF) is shown. 
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Fig. 9. The participant number dependence of multiplicity as obtained from 
the Brueckner G-matrix (GMAT) and a Skyrme force (SMD) are compared with 
the FOPI data [7]. In addition the experimentally derived Harris systematics [45] is 
shown. 
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Fig. 10. The rapidity distributions of Tr"*" as obtained from the Brueckner 
G-matrix (GMAT) and a Skyrme force (SMD) is shown for transverse momenta 
Pt > 0.1 GeV/c. The scaled c.m. rapidity T® = T/Tproj is used. 
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Fig. 11. The tt~ pt spectrum obtained at midrapidity in a An on An reaction at 1 
GeV/A under minimum bias condition is compared to the FOPI data of Ref. [7]. 
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Fig. 12. Same as Fig. 11, but for tt'*'. 
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Fig. 13. The (Te{t^~)/( yE{'^~^) ratio calculated for the same reaction as in 
Fig. 11. at midrapidity under minimum bias condition is shown as a func¬ 
tion of the center-of-mass kinetic energy Tc.m.- The corresponding cross sections 
aE = (Pcr/{pEdEd^l) are obtianed with an angular cut 85^ < ©c.m. < 135°. The 
data points are taken from the FOPI collaboration [7]. The horizontal solid line 
indicates that the observed ratio reaches a constant value of about 0.93 as 

pion energy increases. 







